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Abstract 

We consider a model of bilayer lipid membrane with interdigitation, in which the lipid tails 
of the opposite monolayers interpenetrate. The interdigitation is modeled by linking tails of the 
hydrophobic chains in the opposite monolayers within bilayer as a first approximation. A number 
of thermodynamical characteristics are calculated analytically and compared with the ones of a 
regular membrane without interdigitation. Striking difference between lateral pressure profiles at 
the layers interface for linked and regular bilayer models is found. In the linked case, the lateral 
pressure mid-plane peak disappears, while the free energy per chain increases. Within our model 
we found that in case of elongation of the chains inside a nucleus of e.g. liquid-condensed phase, 
homogeneous interdigitation would be more costly for the membrane's free energy than energy 
of the hydrophobic mismatch between the elongated chains and the liquid-expanded surrounding. 
Nonetheless, an inhomogeneous interdigitation along the nucleous boundary may occur inside a 
"belt" of a width that varies approximately with the hydrophobic mismatch amplitude. 
Key words: bilayer lipid membrane, interdigitation, lateral pressure profile, hydrophobic mismatch 



Electronic address: sergeimoscow@online.ru 



I. INTRODUCTION 



Studying mechanisms of changes in the structure of cell membrane under the adsorption 
of small amphiphilic molecules (alcohol, anesthetics, etc.) is of fundamental interest, as well 
as is important for understanding of the functioning of the cell membranes and embedded 
proteins 111]. One of the drastic changes of the structure is membrane transition into the 
interdigitated phase [2J. While in a regular membrane the thickness of the hydrophobic part 
of a bilayer is approximately twice the length of the hydrophobic tails of the phospholipid, 
an interdigitation may reduce the hydrophobic thickness to the sum of the length of the 
lipid tail and of the small amphiphilic molecule. Understanding possible consequences of 
the interdigitation for the lipid membrane properties is important for e.g. prediction of the 
effects of anesthetics on the functioning of the ion channels embedded in the membrane 

B,Q. 

As a first approximation to the lipid bilayer membrane with interdigitation, in which 
lipid tails from the opposite monolayers interpenetrate, we consider a model with pairwise 
linked tails of the lipids belonging to the opposite monolayers within a single bilayer, Fig. [TJ 
Our model does not allow for a lateral area dilation of the membrane, that may follow the 
interdigitation [lj, but it bears an important property of the interdigitated bilayer in the 
form of constrained meandering freedom of the chains ends in the vicinity of the monolayers 
interface. We found important consequences of this constriction: the entropy of the bilayer 
decreases, the free energy increases, and the lateral pressure profile lit (2) changes drastically. 

The two distinct Ht{ z ) curves, see Fig. [2J can be understood by comparing orienta- 
tional fluctuations of the hydrocarbon segments of the semi-flexible lipid chains in the 
interdigitated- and noninterdigitated lipid bilayers. These fluctuations can be character- 
ized by an orientational order parameter S(z), see Fig. [3] (z is coordinate measuring depth 
inside bilayer), calculated using our model. The fluctuations reach their maximum at the 
monolayers interface inside a noninterdigitated bilayer, because the chains ends are free 
there. Hence, the order parameter drops at z = L, see the dashed curve in Fig. [31 Simul- 
taneously, a maximum of the entropic lateral pressure occurs at z — L, dashed curve in 
Fig. [2j Distinctly, in the interdigitated bilayer fluctuations are significantly suppressed at 
the monolayers interface due to restriction of orientational freedom of the central segments 
by their peripheral neighbors. Hence, S(z) does not drop at z = L, see solid curve in Fig. [3l 
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As a consequence, there is no maximum at z = L in the Ht{z) dependence drawn with solid 
line in Fig. [2] for interdigitated bilayer case. 

We present our analytical results derived in closed form for thermodynamical properties 
of a membrane with linked chains in the weak interdigitation limit: i.e. thickness of the 
hydrophobic part of a bilayer is comparable with twice the length of a hydrophobic chain. 
In this limit we use more complete version of the energy functional entering the membrane 
partition function than developed earlier [sj]: besides the bending energy of a chain confor- 
mation, we included kinetic energy of the lipid chain. We prove that this makes the path 
integral representation of the free energy of the chains uniquely normalizable. 

The plan of the article is as follows. In Sec. II we introduce a microscopic model of a 
membrane with interdigitation and calculate the membrane free energy using path-integral 
summation over the chains conformations. The inter-chain entropic interactions are treated 
in the mean-field approximation. Several thermodynamic moduli characterizing the interdig- 
itated bilayer are derived as well and compared with noninterdigitated case. An increment 
of the free energy per chain due to interdigitation is calculated. In Sec. Ill we calculate an- 
alytically the lateral pressure distribution (profile) across the hydrophobic core of the lipid 
bilayer and make comparison between the cases with- and without interdigitation. Also 
calculated is chain order parameter that characterizes correlations between the orientations 
of the chain segments and clearly demonstrates an increase of the orientation order in the 
interdigitated CclSG clS compared with noninterdigitated bilayer. In the Appendix we dis- 
cuss possibility of inhomogeneous interdigitation along the boundary of the transmembrane 
liquid-condensed domain (raft) embedded in the liquid expanded surrounding in the bilayer 
membrane. We evaluate energetically favorable configuration of such a raft allowing for the 
trade between hydrophobic mismatch and interdigitation-induced free energy increase. 

II. MICROSCOPIC MODEL OF INTERDIGITATED BILAYER 

The interdigitated lipid bilayer membrane is modeled by linking pairwise the tails of the 
chains belonging to the opposite monolayers. Hence, a couple of linked chains is substituted 
by a single semi-flexible string of length « 2L, where L is the monolayer thickness, see Fig. HI 
Correspondingly, conformations of the string as a "trans-membrane" object obey combined 
boundary conditions at the opposite head group regions of bilayer with coordinates z = 
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and z = 2L respectively, as is described below in detail. 

With bending (flexural) rigidity Kf, and with the mean-field approximation accounting 
for entropic repulsion between neighboring couples of pairwise linked chains (see Fig. H]), the 
energy functional of a single string, E t , has the form: 




2 V dz 2 

Here harmonic potential U e ff = BR 2 /2, with self-consistently defined rigidity B, describes 
entropic repulsion between the strings, z is coordinate along the string axis, and R(z) is 
vector in the {x, y} plane characterizing deviation of the string from the straight line, R 2 = 
R 2 + R 2 . The choice of harmonic potential is justified since we assume finite softness of the 
effective "cage" created by the neighboring lipid chains in the limit of small chain deviations. 
A harmonic potential was considered in earlier work for a semi-flexible polymer confined 
along its axis. The first term in Eq. ([I]) represents kinetic energy of the string, p is linear 
density of mass: p = m(CH 2 )N/L, where m(CH2) is a hydrocarbon group mass, N is the 
number of hydrocarbon groups per chain (for numerical estimates we took N = 9, see {(]]). 

The bending energy term in Eq. ([1]) represents the energy of the chain trans or gauche 
conformations. It contains the second derivative over the z coordinate rather than over the 
contour length of the chain. This approximation is valid provided that deviations from the 
z axis are small with respect to the chain length L: 



This limit is opposite to the one considered in the long polymer theory [12|, where the second 
derivative in Eq. (CQ) is substituted by the first derivative in the flexible chain approximation. 

Using the functional Eq. (CD) the chain partition function is found as a path integral over 
all string conformations: 

^((R^R^)) 1 
— ~ I DR x DR * DR v DR v 

(3) 

E[(R x (z),R x (z)) 
exp ^ — J - DR X DR X | = Z 2 r 



k B T 
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The second equality in Eq. ([3]) holds when the membrane is laterally isotropic and the x and 
y deviations can be considered independently. 

To calculate the path integral Eq. ([3]) we rewrite the energy functional Eq. (jTJ) using the 
self-adjoint operator H: 

Et = 2 / ^ + Ri(z)HRi(z)dz) , (4) 

i=x,y 

The operator H is obtained after integrating by parts the expression Eq. (JTJ) under the 
following boundary conditions for the string that models two linked chains belonging to the 
opposite monolayers (the z-coordinate spans from one head group at z = to another at 
z = 2L). The chain angle is fixed in the head group region: 

R'(0) = 0; R'(2L) = (6) 

No total force is applied upon chain at the head group: 

R'"(0) = 0; R"'{2L) = (7) 

These boundary conditions, as well as the energy functional in Eq. (pQ) differ from the ones 
used to describe a single monolayer of a noninterdigitated lipid bilayer (compare jlj]): 

pR 2 {z) K f fd 2 R(z)\ 2 B 



where E™ is the energy functional of a single monolayer, and the motions of the chains in the 
opposite monolayers forming a noniterdigitated bilayer are independent. The total energy 
of a bilayer in this approximation is then twice the energy of a single monolayer: 2 x E™. 
We impose the following boundary conditions for a monolayer: the chain angle is fixed, and 
no total force is applied to the chain at the head group: 

R'(0) = 0; R'"(0) = (9) 

No total force and no torque is applied at the free chain end (i.e. at the monolayers interface 
inside the bilayer): 

R"'(L) = 0; R"(L) = (10) 



Finally, in both cases, the free energy of a bilayer equals F = —kBTln(Z), where Z 
is partition function of a bilayer. Using expressions Eq. (JI]) or Eq. (jSJ) for interdigitated 
or noninterdigitated bilayer respectively, we differentiate the free energy and obtain the 
self-consistency equation in the form (showed for interdigitation case): 

As before Q], we take into account that hydrocarbon chains of lipid molecules are bulky 
objects that possess finite thickness and introduce an incompressible area of the chain cross 
section A (see Fig. [5]). The area occupied by a lipid chain in the bilayer is related to the 
string mean square deviation (R 2 ) by the following formula j^]: 
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cL4 = 7r(R 2 )= (VA-v^To , (12) 



where 5 A is the area swept by the string formed with the centers of the chain cross sections. 
In the text below we imply by chain deviations those of a string described by the R vector. 
The self-consistency equation Eq. fill I) combined with formuli Eq. (1121) permits us to find 
the A dependence of the coefficient of entropic repulsion B and finally derive the membrane 
equation of state in a form of pressure-area isotherm (see Appendix for details). 

To make numerical estimates based on our model of a lipid bilayer we use the following 
parameters values: chain length L = 15 A, chain incompressible area Aq_= 20 A 2 , T = 300 A' 
as reference temperature. The chain flexural rigidity is defined as 9] Kf = EI, where 



E ~ O.QGPa is the chain Youngs modulus 



111 ] and / = Aq/4tt is the (geometric) moment of 



inertia. The flexural rigidity can also be evaluated from polymer theory 12j Kf = kBTl 



"p ■ 

where l p ~ L/3 is the chain persistence length ll| and is the Boltzmann constant. Both 
estimates give approximately Kf ~ kTL/3 at chosen L and at T = T . 

III. INTERDIGITATED BILAYER: THE FREE ENERGY INCREMENT 

The eigenvalues and eigenf unctions of the operator H defined in Eq. (|5]) obey the following 
equation: 

HR n = Kf + BR n = E n R n (13) 
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Solving this equation with the boundary conditions Eq. ©-([Zj) one obtains: 



E n = B + -^,k n = irn/2,n>l; E = B, (14) 



Rn(z) = c n cos(k n z/L),n > 1; R (z) = J —, (15) 



where c n = \/l/L and A n = 2-nL/k n is the wavelength. Several eigenf unctions are shown in 
Fig. El 

Then an arbitrary string conformation, described with the deviation from the z-axis, 
R x (z, t), as well as its energy are expanded over eigenfunctions R n and eigenvalues E n found 
from Eq. ([13]) : 



R x (z,t) = Y,C n (t)R n\ z )'i 



n=0 



(16) 

Rx = C n R n ; E t = - pC n + C n E n 
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n=0 n=0 



The bilayer partition function is then found as the integral over the coefficients of expan- 
sion C n and conjugated momenta p n = pC n in Eq. ( |T6l) 



7 Pi C n E n\ dp n dC n 

Z x = I exp 



n=0 



2pk B T 2k B T J 2nh 

7 rji I 7 rp (17) 



h y E n -"-A ^ 

n=0 ' n=0 



where uj n = \J E r J p. It is important that the latter expression for uj n in the limit o 
free string, B = 0, leads to the well known bending waves spectrum of Euler beam 



a 
fl: 

LU n = \J Elk^l p (with k n = k n /L), as it follows from Eq. (fl4j) and expression for the 
bending rigidity i^/- = £7 mentioned above. Hence, by including kinetic energy of the chain 
into the energy functional E t we obtain correct dimensionless expression for partition sum 
in the Eq. (117)1 . 

Using Eq. (1171) . Eq. ([3D and F = — k B T In Z we find the following expression for the free 
energy of a bilayer with interdigitation in our model: 

n=0 " 



Using then relation Si nt = 



—dFint/dT we find the following expression for the entropy Si nt : 



V hw n J ' ~ [ u n V dT J v T J 
Both expressions are valid provided the motion of the lipid chains at room tempera- 
ture T is classical (not quantum), i.e. : ksT '/hw nmax » 1 and upper cutoff n max in the 
sums is defined by condition that the shortest half-wavelength 0.5X nmax = ^L/k nmax of the 
eigenfunction R nmax is not shorter than the Cif 2- m onomer length of the "chain segment". 

The free energy and entropy of the noninterdigitated bilayer, F non and S non respectively, 
are obtained using the same relations as in Eqs. (fT8l) and (fT9l . but with the corresponding 
change of the frequencies spectrum u n that results from the noniterdigitated bilayer con- 
ditions expressed in Eqs. (lU])- ffTU|) . Using the above relations we calculated interdigitation 
free energy and entropy "cost" as the differences of the respective bilayer free energies and 
entropies in the interdigitated and noninterdigitated cases. Our results are represented in 
Fig. [7J and Fig. [HJ In Fig. [7J the free energy increment (per chain) of the order of 5k bT in 
the intedigitated bilayer with respect to the non-interdigitated one is caused by the corre- 
sponding decrease of the entropy ~ 5/c# (per chain), see Fig. [HJ Location of the entropy 
decrease in the interdigitated bilayer can be found by exploring the chain's orientational 
order parameter S(z) defined as: 

S(z) = 1 -(3(cos*6(z))-l) J (20) 

where 6(z) gives distribution of the tangent angle of the chain across the bilayer. Straight 
(ordered) chain possesses 8 = and S(z) = 1. In the limit of small deviations from the 
straight line 9 < 1 considered in our model the order parameter can be expressed using the 
following relations: 

(cos*9(z)) * 1 - (tg*9(z)) = ((i?'(,)) 2 ) = ^ £ , (21) 

n=0 n 

so that finally we obtain: 

ffW Ml _ ?^ E (22) 



n=0 n 

Calculated order parameter distributions across the bilayer, S(z), in our model with and 
without interdigitation are represented in Fig. [3j The solid line corresponds to linked chains 
(modeling an interdigitation), and dashed line is calculated for non-interdigitated case. It 
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is obvious from the Fig. [3] that main difference occurs at the monolayers interface (z = L) 
inside the bilayer. Free chain ends acquire maximal disorder in this region, while linked tails 
remain quite ordered. Another manifestation of this mid-bilayer ordering phenomenon will 
be seen in the next section in the calculated behavior of the lateral pressure profile inside 
bilayer. 



IV. LATERAL PRESSURE PROFILE AND PRESSURE- AREA ISOTHERMS 
FOR INTERDIGITATED BILAYER 

The equation of state of the lipid chains in the bilayer can be derived as follows: 

!),• 

where P t is the total lateral pressure (or tension), produced by linked hydrocarbon chains. 
Substituting expression for the free energy from Eq. (TT81) into Eq. (123]) one finds: 



OA E v 

n=0 



We may consider P t as an integral of the lateral pressure distribution (profile) function, 
Ht{z), over the hydrophobic thickness of the bilayer: 



P t = J U t (z)dz. (25) 

In order to find out Ht{z) defined this way, it is possible to use the following formal trick. 
Namely, the dependence on aria A of E n arises via dependence of the "potential" B(A), 
that enters operator H in Eq. ([TBI . One may in addition formally consider B(A) as being 
2-dependent function. Then, a well known relation from the perturbation theory [ji]] leads 
to the following equation: 

3E n _ [ 5E n 3B {z)dz _ ( K{z f_m dz< (26) 



OA J 6B(z) dA 1 J OA 
where 1 means unit length. Now, substituting Eq. (126]) into Eq. (1241) we find analytical 

expression for the lateral pressure profile from the relation: 

P. = - / k B Tj2 R f^ 9B f dz = I U t(z )dz. (27) 

J n=0 n J 
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Hence, finally: 



n=0 71 

Calculated in our model lateral pressure profiles for the bilayer with and without inter- 
digitation are presented in Fig. (J2J). It is remarkable, that lateral pressure peak at the 
non-interdigitated monolayers interface, as seen in the dashed curve, disappears in the in- 
terdigitated (linked chains)case. Hence, entropic repulsion between the lipid chains is indeed 
weaker in the region where the entropy related with the chain orientation order is smaller 
(compare with Fig. [3]). 

Next, it is straightforward to check that due to ortho normality of the eigenf unctions 
R n {z) the integral of Ht{z) over dz across the bilayer thickness leads again to the expression 
in Eq. 024p for the total lateral tension P t : 

dB(A)^ 1 , ^^dEn 1 



£^-^£^4, (29) 



dA ^ E n " ^ OA E n 

n=0 n n=0 11 



Pt = ~k B T 

where we used relation that follows from Eq. ([Tj 

In Figure the calculated pressure-area isotherms for interdigitated (linked chains, solid 
line) and non-interdigitated (dashed line) bilayer are presented. It is obvious from the figure 
that lateral entropic repulsion responsible for the lateral pressure in the hydrophobic part 
of the bilayer is weaker in the interdigitated bilayer comparatively with non-interdigitated 
bilayer at the one and the same area per lipid chain and other parameters fixed. 
Differentiation of Pt{A) gives the area compressibility modulus 

K. - -A d £ (31) 

as a function of the area per chain and temperature. The equilibrium condition is found 
by equating the pressure produced by linked chains to the effective lateral pressure in the 
bilayer: 

P t {A{T)) = P eff = j + P HG + P vdW , (32) 



10 



where 7 is the surface tension at the hydrophobic-hydrophilic interface; Pug is the head 
group repulsion of electrostatic origin; P v dw is the pressure arising from the vander Waals 
interactions between chains, etc. We choose P e ff > 7 ~ 70 dyn/cm because attractive 
dispersion interactions between hydrocarbon chains are included in the effective surface 
tension [jjj. At room temperature for a typical lipid bilayer with effective surface tension 
one has: 50 < P e // < 150 dyn/cm jis), Qj- Analytical solution for the total presure in case 
of linked chains (interdigitation): 

^ = 3 -!)'/» ' (2 " 2/S( ^ " 1,2/3 + l) <33) 

Analytical solution for the total pressure in case of not linked chains (no interdigitation): 

Oh T 

It follows from the analysis of these expressions that interdigitation effect on the total lateral 
pressure at a given area is more pronounced at larger areas per lipid (lower pressures) region, 
that corresponds to hihger orientational disorder of the chains. 

In Fig. [10] the temperature dependence of the area per chain in the bilayer is shown. This 
curve is increasing with temperature due to a more frequent collisions of chains. Fig. [TT] 
displays Pt{A) dependence. This curve is decreasing with the area per chain due to the 
following reason: when the chains occupy more space they collide less frequently and produce 
less entropic pressure. 

Now we can verify the exploited approximation of the small chain deviations in the 
bilayer Eq. (T5]). We calculate the thermodynamic average of the chain fluctuation amplitude 
(R 2 (z)) using the relation (B%. y (z)} = J2 n (^n) ^n( z ) an d averaging over C n : 

(R>(z))=k B Tj2^- (35) 

n 

It is worth mentioning that integration of both sides of Eq. fl35|) over z from to 2L 
provides the self-consistency equation Eq. (fTTj) . Since E n oc n 4 , the sum in Eq. (1351) converges 
fast and allowing for the relation R^(z) ~ l/£, we can estimate it as ^2 n l/E n oc 1/B. 
According to Eq. ([TJ), the increase of potential energy associated with the increase of area 
swept by the string from to 8 A is of order B2L5A (the string is formed by the centers of 
the chain cross-sections). On the other hand, it is equal to the work against the pressure 
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P e ff needed to increase the area per couple of linked chains in the bilayer from Aq to 

A: B2L5A m P e ff(A — Aq). From the last equality and relation Eq. (|12p it follows that 
B > P e ff/L. Then we evaluate: 

£^<^-, (36) 

n & n JT e ff 

and find a rough estimate for the upper limit of the small parameter: 

V^R)/2L < (k B T/L 2 P eff )V* = 0.16. (37) 



Finally, we compare the amplitudes of linked chains fluctuations in the bilayer Eq. f[35|) 
and in empty space. For a free couple of chains with flexural rigidity Kf the characteristic 
deviation R° can be evaluated by equating the chain bending energy to k^T . This yields 

R o x ks! L s „ L 2 (3g) 
Kf 



Allowing for Eq. <$S) and Eq. (J3HD , we find ^(R 2 ) / Rq ~ 0.1. 
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APPENDIX A: SOLUTION OF SELF-CONSISTENCY EQUATION 

Here we present the solution of the self-consistency Eq. ( TTTj) and find the analytical 
temperature and area per couple of chains dependence of the lateral pressure produced by 
the linked chains using the equation of state Eq. (|23|) . It is convenient to perform the 
derivations in dimensionless parameters, 

a = A/A ,b = ^-B, (Al) 
and to introduce the auxiliary parameters 

k n = {itn/2)\n >l;v= f^" (A2) 

where L ~ lb A is the chain length, A ~ A 2 is the "incompressible area" of the chain 
cross section, and the chain flexural rigidity Kf = kBTL/3 at T « Tq = 300 K. Using these 
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estimates we obtain v = 0.009. 

In the introduced notations Eq. (lAlj) and Eq. (1A2I) with E n defined in Eq. (1141) the 
self-consistency equation Eq. (TiTj) acquires the form 

i^^r = ^-lf. (A3) 



The terms in the sum on the left hand side of Eq. (I A3[) decrease fast with growing n and 
we can use integration instead of summation over n. For example, tt the effective tension 
P e ff = 70 dyn/cm, we have b m 10 3 , while kil 4 ~ 9 • 10 4 . In this regime we can solve 
Eq. (1A3|) analytically by substituting summation over n with integration, which yields 

\- 1 j 1 r dn 1 

2^ A _|_ M ~ O / J, -I- fc 4 .. , AJ, 3/4 I J 



x b + K 2 7^ 6 + ki +1 ^263/4 
where we took integral using complex functions theory, and k n is defined in Eq. (1A2|) . 
In case of membrane with no interdigitation (see Fig. Hj) Eq. ( 1A3I) takes the form: 



^ + E^II = ^-1) 2 - (A5) 



b ^b + k* 

n=l " 

Since b ~ 10 3 (-P e // = 70 dyn/cm), we (approximately) integrate over n, so that in case of 
membrane with no interdigitation Eq. ( 1A4j) takes the form: 

v 1 „ 1 r ^ _ 1 fA61 

^ & + & 4 2j_ 00 b + k* +1 2V2b^' 1 J 

where fc n is defined as in 5j: /c n = 7rn — 7r/4. 

In both, Eq. (1A4|) and Eq. (1A6I) . we omit 1/6 as 6 ~ 10 3 , which leads to the same b(a) 
dependences in both cases: 

(A7) 



4w 4 / 3 (v^- l) 8 / 3 ' 

which is then used in the equation of state Eq. ( 1231) . As a result we find the expression 
for the lateral pressure produced by the linked hydrocarbon chains Eq. (123)) 
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F L p, k B T 


i 


0.61 


2 


2.43 


3 


5.46 


4 


9.71 


5 


15.17 


6 


21.84 



TABLE I: Energy of hydrophobic mismatch {Flp see Eq. IBlj) as a function of raft thichness (L is 
a lipid membrane monolayer thickness, i.e. constant) 

APPENDIX B: TRANS-MEMBRANE STRUCTURE OF RAFTS 

Here we make some estimates to find out if there could be free energy gain related with 
interdigitation along the perimeter of hydrophobically mismatched regions in the membrane, 
e.g. for the case of a cluster of the ordered lipids surrounded by the liquid membrane 
"sea" (model of a raft). We evaluate energy of hydrophobic mismatch using the following 
expression: 



LP 



77 2t" 

-(K d rl + K a ){^ + l){d p -Lf 
4 r 



(Bl) 



where r p ~ 25A [16| is a protein radius, rnjy 10A (13] is a characteristic scale of deforma- 
tion (see Fig. [12]) . ~ 15 • 10 1 4 erg /cm 4 18] is a dilation modulus, K g ~ 35 erg /cm 2 [3] 
is the modulus, characterizing the energy cost of producing a gradient of bilayer thickness 
(it includes the energy of increasing the area of chain- water interface), L is the (equilibrium) 
monolayer thickness, d p is the protein monolayer thickness. One can see Flp values with 
given values of parameters and various d p — L, i.e. raft-bilayer thickness micmatch. 

Thus, we found that Flp ~ lOkpT whereas energy cost of interdigitation in the ro 
area (see Fig. [T2l) is Fi nt ~ 100k pT (one can easily obtain this value by counting chain 



quantity in the r area and multiply its by AF 



perchain\ 
int ) 



This means, that, within our model, 



interdigitation itself does not help to decrease the hydrophobic mismatch energy of the 
cluster (raft). 
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FIG. 1: Model of membrane with interdigitation we based on. 




FIG. 2: Lateral pressure distribution in the hydrophobic core of the bilayer. z is coordinate along 
the chain axis normalized by the monolayer thickness L and spanning from one head group {z = 0) 
to another (z = 2L). The parameters for the lipid bilayer are as follows: monolayer thickness 
L = 15^4, area per chain Aq = 20A, chain flexural rigigity K ~ kTL/3, temperature T = 300-ftT, 
total monolayer pressure P = 70 dyn/cm. 
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FIG. 3: Order parameter in cases of linked chains (see Fig. [TJ) and no interdigitaion (see Fig. [4]). 



-t 

single chain 
+ U eff 

bilayer 

FIG. 4: Model of lipid membrane in the mean-field approximation: we substitude interaction 
between neightboring chains by an effective quadratic potential. 




FIG. 5: Hydrocarbon chain as a flexible string of finite thickness. R(.z) is the vector characterizing 
the deviation of the center of the chain cross section from the z axis, |R(z)| = \J R x { z ) + Ry{z)\ A) 
is the incompressible area of the chain cross section; A = ir (R- 2 ) is the area swept by the centers 
of chain cross sections; A is the average area per lipid chain in the bilayer. 
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FIG. 6: The eigenfunctions R(z) of the self-adjoint operator H for the boundary conditions Eq. © 
and Eq. (|7|). Other parametrs are as in Fig. [2l 
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FIG. 7: Not normalized free energy (per chain) difference of membranes with linked chains and 
with no interdigitation. Linked chains membrane "cost" more free energy. 
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FIG. 8: Entropy (per chain) difference between membranes with linked chains and with no inter- 
digitation. Entropy of linked chains membrane is lower that of no interdigitation membrane. 
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FIG. 9: Total lateral pressure comparison. With the same are per chain, membrane with linked 
chains produce less pressure than no interdigitation membrane. 

2,48-, 
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FIG. 11: Calculated total lateral pressure Pt produced by linked hydrocarbon chains as a function of 
area per chain at two temperatures T\ (solid line)< T2 (dashed line). Lateral pressure is normalized 
by hsT/A®, area is normalized by Aq. Other parameters are as in Fig. ([2]). 
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